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By implementing an additional heavy quark—antiquark pair production trigger in a multiphase transport 
(AMPT) model, we study the effect on anisotropy flows of identified particles with a focus on charged particles 
and quarkonium (J/W and Y). A systematic increase in the collision rate for active partons in the AMPT model 
with such an implementation has been observed. It leads to a slight increase of identified particles anisotropy 
flows as a function of transverse momentum (pr) and rapidity, and gives a better description of the experimental 
data of elliptic flow toward larger pr. Our approach provides an efficient way to study the heavy quark dynamics 


in the AMPT model at LHC energies. 
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I. INTRODUCTION 


In nuclear—nucleus collisions at the Relativistic Heavy Ion 
Collider (RHIC) and Large Hadron Collider (LHC), a sys- 
tem of color deconfinement, known as Quark-Gluon Plasma 
(QGP), is expected to be produced [1, 2]. QGP is an al- 
most perfect fluid in the sense that its shear viscosity to en- 
tropy density ratio is close to the minimum value [3, 4]. Be- 
cause the volume and lifetime of the QGP may be much 
larger and longer than those given by the confinement scale 
(~ 1/Agcp), it is expected that the collective flow of par- 
tons should arise and have observable consequences that pro- 
vide information on the dynamics in the QGP [5]. Studies on 
the anisotropic flow of particle production have contributed 
significantly to the characterization of the system created at 
RHIC and LHC [5]. Using a general Fourier series decompo- 
sition of the azimuthal distribution of the emitted particles, 


dN =~ 
Ip 1+2 X vn cos[n(y — U,)], d) 


n=1 


anisotropic flow is quantified using coefficients v, and the 
corresponding symmetry planes Y, [6], where y represents 
the azimuthal angle of the emitted particles in the momentum 
space. Systematic measurements of the correlations and fluc- 
tuations of v, coefficients and Y, have been performed in 
Refs. [7—13]. Comprehensive comparisons with theoretical 
model calculations provide critical information on the event 
average initial density distribution in the nuclear overlap re- 
gion, as well as its connection to QGP dynamics [14-19]. 
Charm and bottom quarks are important probes for the 
QGP. They were mainly produced in the initial state of the 
hard-scattering process prior to QGP formation and survived 
throughout the evolution stage of the QGP. Quarkonia are 
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bound states of heavy flavor quark—antiquark pairs, which are 
believed to be suppressed because of the color screening ef- 
fect in the deconfined interior of the interaction region [20]. 
They offer a complementary way to study the interaction of 
the heavy flavor with the medium and, thus, independently 
shed light on the properties of QGP [21-23]. Dozens of mea- 
surements in this direction have been carried out, and rich 
physical information has been delivered [24-31], with a clear 
future direction [21]. 

To study the characteristics of the QGP and its evolution to- 
ward chemical and thermal equilibrium, it is necessary to de- 
rive transport equations for quarks and gluons when they are 
in equilibrium and pre-equilibrium [32]. There have been sev- 
eral studies on deriving relativistic abelian and non-abelian 
transport and constraint equations for the relativistic gauge 
covariant Wigner operator for color particles interacting via 
SU(N) gauge fields. This indicates that viscous hydrody- 
namic models can describe the transverse momentum spectra 
of light-flavor hadrons and distribution in the azimuthal angle 
for low pr [16, 33]. On the other hand, transport models have 
been applied to study the evolution of QGP medium [34]. 
Owing to its relatively small cross section, a large event sam- 
ple is required to study charm or even bottom quark transport 
in the QGP within hydrodynamics calculations or transport 
model simulations. In this study, we implement an additional 
heavy quark—antiquark pair production trigger in a multiphase 
transport (AMPT) model to study the dynamics of charm and 
bottom quarks in a dense light quark soup. This significantly 
improves the charm or bottom quark simulation efficiency, 
and it will be closer to the data at LHC energies because of a 
large charm quark production cross section [21]. 


Il. TRANSPORT MODELS OF HEAVY FLAVOR 


We employed the AMPT model with the string melting 
version (v2.265) to simulate the collective motion of heavy 
quarkonium in Pb—Pb collisions at LHC energies. The AMPT 
model is a hybrid model consisting of fluctuating initial con- 
ditions, parton elastic scattering, quark hadronization, and 
hadronic interaction [34]. It has been used extensively to de- 
scribe the dynamics of heavy-ion collisions at RHIC and LHC 


energies in some recent publications [35—42]. For a detailed 
description of the AMPT model, please refer to Ref. [34] and 
Ref. [43] for a summary of recent developments. 

In the AMPT model, the initial production of heavy quarks 
is handled using the HIJING two-component model [44]. It 
includes pair productions (q + g > Q+Q,g+9 > Q+ 
Q) and gluon splitting (g — Q + Q) with a cross section 
in perturbative quantum chromodynamics (QCD) at leading 
order: 


do 
Zi fal£ı, 
dpždyıdy2 oe talea) 
dg >Q 
x2 fo(w2,uR) X —F (2) 


Here, yı and yo are the rapidities of the two produced partons, 
the K factor aims to account for higher-order corrections of 
heavy quark production, a and b refer to the types of interact- 
ing partons in the initial state, x denotes the nucleon momen- 
tum fraction taken by the interaction parton, up represents the 
factorization scale, f, is the parton distribution function of 
parton type a in a nucleon, and 07> is the cross section 
for parton types a and b to produce the heavy quark pair [44]. 
Because the cross section of heavy quark production, espe- 
cially for the bottom quark, is several orders of magnitude 
smaller than that of light quarks in heavy-ion collisions at 
LHC energies, it calls for a trigger for heavy quark produc- 
tion to increase the simulation efficiency. The algorithm was 
originally developed by Wang and Gyulassy [44] and was ex- 
tended to study the bottom quark dynamics in the QGP in this 
study. We refer to the details in Ref. [44] and Ref. [45] and 
do not repeat them here. 

The final state parton interactions can be described by the 
equations of motion for the quark and gluon Wigner opera- 
tors [46, 47], which are closely connected with the classical 
distribution function. For example, the relativistic gauge co- 
variant Wigner operator for spin-4 particles in coordinate and 
momentum space [46]: 


4 
w (xp) = J a 


ec PUa (x) e39 Ds 


Be 24 Poy (x), (3) 


where D” = 0" + ig A is the covariant derivative, and W is 
a matrix in spinor (4 x 4 components) as well as in the color 
space (N x N components). 

By using the properties of the link operators, one can de- 
rive the quantum transport equation for the QCD Wigner op- 
erator [46, 47]. Under semiclassical and Abelian dominance 
approximations, in the weak-field limit, the equations of mo- 
tion become a set of equations that resemble classical trans- 
port equations [48]. This motivated the study of the final state 
interaction by solving the Boltzmann equations for quarks 
and gluons. If only two-body interactions are considered, the 
Boltzmann equations can be reduced to the following: 


p"O, f(a, p, t) x J oflesps.t)f(ea,pa.t), (4) 


where f(x, p,t) represents the distribution function of the 
parton at time ¢ in the phase space, and ø is the cross sec- 
tion for the two-body scattering of partons. Zhang’s parton 
cascade (ZPC) model was developed to study the final-state 
parton interactions, which solved the Boltzmann equations 
for partons using the cascade method [48]. In this model, 
two partons scatter when their closest distance is less than 
,/o/n. The scattering angle is determined by the differential 
parton—parton scattering cross section. Because the cascade 
prescription for solving the semiclassical transport equation 
does not assume any particular form of the phase space dis- 
tribution, it can be used to study systems far from local ther- 
mal equilibrium, and it automatically takes into account the 
finite mean free path effect [49]. The current study focuses 
on heavy quark dynamics at LHC energies, where the bottom 
quark may be far from thermalization and the charm quark 
may be close to thermalization [27, 28, 31]. 


Ill. RESULTS AND DISCUSSION 


To follow the parton scattering history in the AMPT model, 
we define Nscat. as the number of two-body elastic scatterings 
suffered by a parton. The Ngcat. distribution of the active par- 
tons (Nscat. > 0) is shown in Fig. 1]. The distributions from 
the AMPT model with an enhanced heavy quark trigger are 
flatter than those in the normal AMPT calculation, which sug- 
gests a larger scattering possibility for partons with enhanced 
heavy quark triggers in the event. The effect is rather small 
for the charm quark scenario and more significant for the bot- 
tom quark scenario. 

Figure | also shows the approximate average transverse ra- 
dial velocity (81) = (tf. - P1) of freeze-out partons with 
different flavor in the AMPT model as a function of Ngcat., 
where f, and p, are the transverse radial position and mo- 
mentum unit vectors. The 3, is not exactly zero at Nscat, = 0 
because partons can form only after a finite formation time 
over which the displacement of a parton depends on its mo- 
mentum [34]. One sees that 6, of heavy quarks are larger 
in magnitude and almost independent of Nscat., whereas it is 
smaller in magnitude and drops for light quarks as Nscat. in- 
creases. This is understood because heavy quarks are formed 
earlier, and they are more likely to inherit the direction of 
mother partons, which will give them larger (3, ). Compared 
with light quarks, it is harder to push the heavy quarks away 
from their original trajectory, their (81) are independent of 
Necat.- AS the mid-rapidity hadron density of a Pb—Pb colli- 
sion at ,/Snn = 5.02 TeV reaches as high as a few thousand, 
depending on the collision centrality [50], triggering on an ad- 
ditional heavy quark—antiquark pair production has a negligi- 
ble effect on the hadron density, but may affect the evolution 
dynamics significantly. We use these trigger events to study 
the properties of QGP matter created in Pb—Pb collisions at 
LHC energies. 

Figure 2 presents v2, v3 of a charged hadron in mid- 
central Pb—Pb collisions at \/snn = 5.02 TeV from AMPT 
model. To do a fair comparison with experimental data [51], 
the results are calculated using the two-particle correlation 
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Fig. 1. AMPT simulations of: (a) normalized parton scattering dis- 
tributions of active quarks; (b) mean transverse velocity of freeze- 
out quarks in Pb-Pb collisions at \/snn = 5.02 TeV and 5%-60% 
centrality, where normal, w/ c-trig., and w/ b-trig. represent normal 
AMPT, normal AMPT with triggering an additional charm quark— 
antiquark pair production, and normal AMPT with triggering an ad- 
ditional bottom quark—antiquark pair production, respectively. 


method [52]. The result from the reaction plane method, as 
described in Eq. (1), has been checked and is found to be con- 
sistent with Fig. 2. Overall, the AMPT calculations with the 
normal, enhanced heavy quark trigger scenario reproduce the 
experimental data reasonably well. Larger elliptic flow is seen 
in the enhanced heavy quark trigger scenario. The enhance- 
ment is approximately 5% and pr independent in AMPT with 
bottom quark trigger (cf. the bottom panel of Fig. 2). In- 
cluding the discussion in Fig. 1 where the AMPT model with 
enhanced bottom quark—antiquark pair production presents a 
larger scattering possibility for active quarks, which tends to 
generate larger elliptic flow and is seemingly closer to v2 data 
as pr increases. 

Figure 3 highlights the v2(pr) of quarkonia in middle and 
forward rapidity of Pb—Pb collisions at \/snn = 5.02 TeV. 
The uncertainties in calculations are determined by statistics 
in different centrality and different rapidity bins. Our results 
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Fig. 2. Charged particles elliptic flow (v2) and triangular flow 


(v3) in mid-central Pb-Pb collisions at \/snn = 5.02 TeV from 
the AMPT model with normal, enhanced heavy quark production 
trigger settings. The filled points represent experimental data for 
comparison [51]. 
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Fig. 3. Heavy flavor quarkonia v2(pr) of mid-rapidity and forward 
rapidity in Pb-Pb collisions at ,/snn = 5.02 TeV from the AMPT 
model with enhanced heavy quark—antiquark pair trigger settings. 
The centralities were chosen to be the same as those in the experi- 
mental measurement [31, 53]. 
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Fig. 4. AMPT model simulations on hadron v2(7) in Pb-Pb col- 
lisions at \/snn = 2.76 and 5.02 TeV. The centralities and pr bins 
were chosen to be the same as the experimental data [31, 54]. Dif- 
ferent model configurations were applied in the charge hadron calcu- 
lation, whereas only the heavy quark—antiquark pair trigger settings 
were used for studying J/W and Y. 


are consistent with the experimental data considering current 
uncertainties. Overall, AMPT simulations on v2 of Y are 
smaller than those of J/Y. One should note that in the cur- 
rent version of the AMPT calculation, regeneration of J/W or 
dissociation of Y are missing, only continued two-body elas- 
tic scattering is involved, which seems to generate the finite 
v2 of quarkonia. The magnitude difference between Y and 
J/W may be due to the mass of different quarks associated 
with each calculation. We have learned that it is also one of 
the key parameters in hydrodynamics calculations [16, 33]. 
To study this further, we examined the differential rapid- 
ity (pseudorapidity) dependence on v2, as shown in Fig. 4. 
Because v2(7) data for a charged hadron in Pb-Pb collisions 
at ,/s~n = 5.02 TeV has not yet been published, we chose 
the data [54] at \/syjn = 2.76 TeV for comparison. Thus, 
the AMPT calculation for the charged hadron is performed on 
/8nN = 2.76 TeV accordingly. Our calculations with normal 
or enhanced heavy quark—antiquark pair production scenar- 
ios describe the data fairly well. Although the uncertainties 
presented in the measurement do not allow us to distinguish 
between different scenarios, the difference in the calculations 
is of a few per cent. The value of Y və is positive and at a 
maximum at mid-rapidity, and close to zero at higher rapid- 
ity. The calculation appears to agree with the data at forward 
rapidity [31]. The value of J/W ve is larger than Y and ex- 
hibits similar rapidity dependence. However, it underpredicts 
the data at forward rapidity. This deviation may originate 
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Fig. 5. AMPT model calculations on hadron v: (7) in Pb-Pb colli- 
sions at \/Snn = 5.02 TeV. The charged hadron result is scaled by a 
factor of 10 for comparison with the experimental data [55]. 


from the missing dissociation and regeneration mechanisms 
for heavy quarkonium in the AMPT. Further studies are re- 
quired to implement these effects in the AMPT model. 

The enhanced heavy quark—antiquark pair production event 
can also be applied to study the direct flow (v1) of quarko- 
nium. This represents a new probe for studying the strong 
electromagnetic field effect of hot QCD matter via vı mea- 
surement, particularly on the difference in v;(7) between the 
particle and its antiparticle [56-58]. Measurements of the 
open charm meson vı(ņ) show a difference in the slopes 
of D° and D°, but with a different sign from RHIC to 
LHC [55, 59, 60]. The origin of this physics is still debated. 
Because the electromagnetic field effect has not been imple- 
mented in the AMPT model, we consider the vı (n) of J/V 
and Y instead of the vı difference between the particle and 
antiparticle in our current framework. In this aspect, we used 
the coordinate system of the model with the z-axis along the 
impact parameter direction (b) of each event in the transverse 
plane, with the center of the projective and target nuclei at 
the transverse coordinates (b/2, 0) and (—b/2,0). The direc- 
tion of the projectile momentum specifies the z-axis; conse- 
quently, the total angular momentum of the system is along 
the —y direction. Figure 5 shows our results for J/WV and 
Y together with the charged hadron in Pb-Pb collisions at 
V/SnN = 5.02 TeV. It can be seen that quarkonia present 
an order of magnitude larger v;(7) than charged hadrons, 
presumably because of the early formation and larger mass. 
In comparison with the data [55], the v1(7) of the charged 
hadron from the AMPT model is on the correct sign but an 
order of magnitude lower. Because the dynamics that gener- 


ates vı is quite complicated, it may depend on the equation 
of state, mean-field potentials, particle scattering, shadowing 
from spectator nucleons, and tilt of the hot dense QCD matter 
in the x-z plane [58]. Our results rely on two-body parton 
elastic scattering, and naive quark hadronization can be un- 
derstood. Future studies with more realistic dynamics are of 
interest. 


IV. SUMMARY 


In summary, we have implemented an additional heavy 
flavor quark—antiquark pair production trigger in the AMPT 
model to study their dynamics in hot and dense media created 
in high-energy heavy-ion collisions. In this study, we have 
focused on the heavy flavor quarkonium anisotropy flows in 
Pb-Pb collisions at LHC energies. The new implementation 
presents a systematic increase in the flow coefficients as a 
function of pr and 7 in comparison with the result from the 
normal AMPT model and provides an efficient way to study 
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